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ABSTRACT 


A new technique for observer design which requires the 
Scion ota minimax problem is presented. An algorithm 
for solving the minimax problem is proposed and applied to 


the design of full-order and reduced-order observers. 
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I. INTRODUCTION 


risco scono vistems, the input vector to a plant 
we cena l unction of the current plant state vector. This 
means that one requires knowledge of the plant state vector 
dine current time to accomplish system control. In most 
Rien, however, the entire state vector is not available 
for direct measurement. Hence, it is desirable to recon- 
eve, or estimate, the unknown plant state vector. The 
device which performs this reconstruction of the plant state 
vector is called an "observer." In this thesis, after 
presenting necessary background material on the theory of 
observers, which is primarily the work of Luenberger [4], a 
new approach to observer design is systematically developed. 
This new method uses a quadratic performance measure which 
Bro be minimized with respect to the observer gain matrix 
G. 

samp lic, Only linear time-invariant continuous- 
time systems are considered throughout this thesis. 

Maca licor nmproposed as suitable Tor digital computa- 
mei Using appropriate Subroutines and is applicable to 
Bacher-order systems. In order to illustrate the new 


Pcie. numer Cal examples are presented. 








CESOU E O ERVER CHEORY 


A. THE NEED FOR OBSERVERS 

To design a feedback control system it is often desired 
to Know the entire plant state vector. But, in most cases, 
WRP entire state vector cannot be measured, hence a suitable 
estimate to the plant state vector must be determined. This 
creates the need for a system called an "observer" which 
produces an estimate of the Piave tor, 

Recunine char 1b as possible to Build a system which can 
estimate some constant linear transformation of the system 
Rate vector, then, provided that the transformation is 
miertible, it is possible to reconstruct the state vector 
itself. This is the fundamental idea upon which observer 


theory is based. 


E FULL-ORDER OBSERVERS 
First consider the construction of a full-order observer 


moma iree system-described by 


x(t) 


li 
[> 


x(t) (1) 


NO) 


li 
| 


x(t) (2) 


Wire X(t) is the nxi state vector, y(t) is the mxl output 
vector, À is an nxn constant matrix and C is an mxn constant 
matrix, and n > m. Thus, an m-dimensional output vector y(t) 
is available from measurements, but we want to know the 


cure spate vector xt) of the plant. Since x(t) cannot be 





measured directly, one must obtain an estimate x(t), that is 
a good approximation to the original plant state x(t) for 
eel times In what follows it is assumed that the plant 
Sale ribed by (1) and (2) is observable. | 

One way to estimate the state vector x(t) might be to 
construct a system of the same dynamic order and structure 


as the plant, given by 
EN ao) (195) 


amar ive it in parallel with the plant. If this is done 
Ehe System state, z(t), will be a linear transformation of 
Wi mplant state x(t). If the initial states z(0) and x(0) 
are equal, then z(t) and x(t) remain equal for all t 2 0. 
Pevarly, however, z(0) is not equal to x(0) and the differ- 
ace between the two states a(t) and x(t), can be determined 


By Torming 


z(t) =- x(t) = A [z(t) - x(t)] (4) 
ime solution of (4) is 
z(t) - x(t) = em! [2(0) - x(0)] (5) 


mus) the error zit) — x(t) depends on the matrix A. If the 
Emma Aos a stability matrix, i.e., all the eigenvalues of 
the matrix A have negative real parts, the error approaches 
ec as cime goes Lo Infinity. 

A better alternative is to construct a linear system 


which contains a model of the plant state equation but is 








eo driven by an output error signal, y(t) = C z(t), that 


19, 
z(t) = A z(t) + G[y(t) - C z(t)] (6) 


or 


Pe Ce Co E E ct) (7) 


This system has a very desirable feature: the observer gain 
meet G can be selected to determine the rate at which the 
observer state z(t) approaches the plant state x(t). If the 


state z(t) becomes equal to the plant state x(t) at some time 


E E ente Correction term, y(t) — C z(t), is zero for 


IL 


t = t. and the state z(t) will remain aligned with the plant 


Ji 
spate x(t) for t > tı. 

The system described by equation (6) or (7) is referred 
to as a "full-order observer" for the plant described by (1) 
and (2). The gain matrix G can be selected by the designer. 
RENO select the gain matrix G is theyproblem of observer 
design. 

The above development can be easily extended to forced 


Aereas by including the Input ult) in the observer equation 


(6) or (7). Consider the system described by 


Il 


x(t) = A x(t) + B u(t) (8) 


y(t) 


Ji 
{© 


x(t) (9) 





Mere u(t) às the rx! input vector and B is an nxr constant 
biisiviriibpution matrix. Tne observer for this plant is similar 
Nai or the free system, but the input vector u(t) must 


be included. The appropriate full-order observer system is 


z(t) = A z(t) + Bult) + Gly(t) - C z(t)] (10) 


or 


z(t) = (A - GC) z(t) +B u(t) +G y(t) (11) 


me configuration of this full-order observer is shown in 


Hic. l. 





FIGURE 1. FULL-ORDER OBSERVER 








C.  REDUCED-ORDER OBSERVERS 

The full-order observer possesses a certain degree of 
redundancy. The reason for unas tS "that tie Tull- 
Bullen Observer constructs an estimate of the entire plant 
meee. DUG che output of the plant, which is made up of 
linear combinations of the states, is available by direct 
measurement. Hence, it ís actually unnecessary to construct 
a system like the full-order observer to estimate the entire 
plant state. 

ines redundancy can be eliminated by reducing the dynamic 
order of the observer system to (n - m), the difference 
between the dimensions of the state vector x(t) and the 
Rui vector vit). The observer which is reduced to its 
minimal dynamic order (n - m) is called a "reduced-order 
observer." The structure of a reduced-order observer is 


WMiliostcrated in Fig. 2. 
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FIGURE 2. STRUCTURE OF REDUCED~ORDER OBSERVER 


iiMewsuarcineg point for consideration of reduced-order 


SUsSemvyers is apain the plant characterized by 


x(t) x(t) + B u(t) di) 


I 
|> 


11 
O 


y(t) x(t) (13) 


TO 





Eoc output distribution matrix C has dimension m. By intro- 
wanne a change of coordinates, it can be assumed that the 
matrix C takes the form C= | I O nem stage 
vector X(t), is selected sothat its first m components are 
equal to the plant output vector y(t) and the remaining 

(n - m) components, denoted by w(t), can be selected arbi- 
trarily but in such a way that the states y(t) and w(t) 

are linearly independent. Thus, the new state vector X(t) 


Sam De Written in partitioned form as 


and x(t) is related to the original plant state vector x(t) 


Poem nonsingular matrix P, l.e., 


ie) 


I 
has, 
rs 
~^ 
ct 
ns 
Il 
Pa 
^~ 
CT 
n 


The matrix C here is the one given in (13) and D is an 
aoa Crary matrix selected in such a way that E. is non- 
eeagoular., The new state equations have the form 


y(t) y(t) 
ASA Teen = A! _ æ — = A, (14) 


w(t) w(t) 


inesmaerices hl and BY in (14) are not the same as A and B 


in (12). These two sets of matrices are related to one 


TL 








another through P. Equation (14) can be rewritten in 


parcicioned matrix form as 


y(t) Au, Aap y(t) B, 
n~ | KÎS SS nan nî n] J ~~ < | A = - late) (i5) 
us. ee a Bar 
; | 
Or 
y(6) = Ay," y(6) AW w(t) + By! u(t) (16) 
wet) = Aj! y(6) + Ago! wet) + B,' ult) (17) 


WWw-fivector y(t) is available for measurement, and if it is 
bilWicrentiated, y(t) can also be determined. Since u(t) is 


also available, (16) provides A,.,' w(t) which plays the role 


T2 
of a set of measurement equations for the system (17) which 
has w(t) as íts state vector and Ax)! y(t) + B,' u(t) as its 
NE vector. The next step is to construct a full-order 


observer of order (n - m) for the plant of equation (17) 


Using 


A 


Rea) 


as the measurement. Thus, the observer for the system 


described by (16) and (17) can be expressed as 


2 





A A 
figs) = At GA) lt) + A, "y(t) + Bt u(t) 
+ G [y (€) - Ay,’ u(t) - 6 B,* ute) (18) 
or 
A A 
meee) = (AOD GAS wt (AL G Ay yi) 
+ (B," - G Bi") u(t) + G y(t) (19) 


Dun ve w represents the estimate of w. The gain matrix G 
in (18) or (19) can be selected so that the observer system 


Î ! son 
22 G A,,' has arbitrary eigenvalues 


provided that (A,,', A,,') is observable. It is known [4] 


that if (C, A) is observable, then (A,,', A,,') is also 


de iii cient matrix A 


sbéervable. The configuration of the reduced-order observer 


shown in Fig. 3. 


E 


Ut) , î 
ji | Au — û A 


> 


u(t) 





FIGURE 3. REDUCED-ORDER OBSERVER WITH DERIVATIVE 


o 





The differentiation of y(t) in (19) may be troublesome 
Mommer horn physicatliy; however, it can be avoided by modi- 
fying the observer state equation (19). The resulting 


mođified observer state equation can be written as 


BR = (ro SGA!) alt) + (Aa! - 04,5!) G y(t) 


WA. a I ue) 
(20) 
with 


w ey e) (zus) 


and the corresponding observer diagram is shown in Fig. 4 


which is equivalent to the observer in Fig. 3 at the point 


W. 





u(t) 


B.-GB. | | Án -SAn 


FIGURE 4. MODIFIED REDUCED-ORDER OBSERVER 
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D. THE CONVENTIONAL METHOD OF OBSERVER DESIGN 
It has previously been shown that the full-order observer 


for the linear time-invariant plant 


x(t) = A x(t) + B u(t) 
Cec) 
y(t) = C x(t) 
is characterized by 
z(t) ee Momo cL) + acy (t) +8 B ult) (23) 


If (C, A) is observable, the observer's eigenvalue can be 
M@eaveQ arbDitrariiy by selecting the gain matrix G: The 
problem is how to select the eigenvalue locations for the 
Bbserver. | | 


The observer error e(t) is defined as 
A 
e(t) Ê z(t) - x(t) (24) 


Differentiating equation (24), substituting the expressions 


rêr z(t) and x(t), and solving for e(t) gives 


e(t) EG GC) ein) (259 


ES Cine; 


p £ 


ã& - & € 


then the error equation for the observer can be written as 


e(t) = E elt) (26) 





Tne solution to this equation is 
RE 
e(t) = E * e(0) (27) 


Mins? the dynamic response oí the error between the observer 
Enche plant state is determined by the matrix F = A - GC. 
It is desired that the error e(t) eventually approach zero 
Ml e(0). Thus, If GG is selected so that F = A - G € 

is a stability matrix then the observer state z(t) approaches 
tue plant state x(t) at a rate controlled by G. One approach 
to observer design is to select the desired eigenvalues and 


then determine G to obtain these eigenvalues [4]. 


16 





III. A NEW METHOD FOR OBSERVER DESIGN 


RR PROBLEM DESCRIPTION 

AN alternative procedure for observer design can be 
based on a performance measure for the error system (26) 
Emis cting Oi the integral of a quadratic form in the 


sbiserver error, that is, 


Es 


el(1) Q e(t) dr (28) 


J= f 
to 


Mere 9 is a real symmetric positive definite constant matrix 
and e(t) is a function of the gain matrix G, which is to be 
selected to minimize this performance measure. It is shown 


in the Appendix that equation (28) can also be expressed as 
(G) e (29) 


where A Rear ne Dia error ol Lhe observer, a is the 
transpose of e, and Ve) is the value at t = 0 of the real 
symmetric positive definite solution V(t) of the matrix 


differential eouation 
V(t) = -P Vit) - VIT) F-Q (30) 


Subject to the boundary condition 


O (31) 


Dei 





merernavively, defining t = ve - T, this equation can be 


expressec as 


V(t) =F V(t) + V(t) F+Q (32) 


Tre the initial condition 


NEN (35) 


From equation (29) it is seen that J depends on the 
initial error of the observer, Co: and on the matrix E 
eS implicitly dependent on G through equation (30). 
Finding a solution to equation (30) requires solving 
n(n + 1)/2 differential equations (since V is a symmetric 


Meer aX) . 


B. CONSIDERATION OF THE INITIAL CONDITION 
ER desired to select the observer gain matrix G to 


minimize the performance measure 


Up 


el (rt) Q e(t) dı = 25 (G) 


qe e, 
O 


t 


To calculate the performance measure, a priori knowledge of 
the initial state of the error system is needed. The initial 
error state is not known, but a reasonable approach is to 
take the worst initial state in some bounded region. It will 
be assumed that the initial observer state is 0 and that the 


possible initial error states lie within a hyperspherical 


18 





region with center at the origin and having a specified 
radius p. With this assumption the performance measure can 


be written as 


t 
tom 
ED mess | S e (1) @ e(t) at 
p t ERN 
e 
_ è, 
leoll< e 
= max Te Cc) ment (34) 
= Co: =o 
Ez 
lle ll < e 
where || || denotes the Euclidean norm. It is easily shown 
that the maximum occurs where le || = p , moreover, without 
lees Ol generality it can be assumed that p = 1. That is, 


the maximum can be considered to occur on the unit hyper- 
sphere. This is because, for any e, such that |le || = 1 and 


aayeracdius p > 0 , 


(35) 


li 
ES 
es 


i.e., the maximum performance measure on the hypersphere of 
: C 2 
radius p is p times the maximum performance measure on the 


unit hypersphere. Hence, it is desired to minimize 


19 








a max J e’(r) ar eim) dr 

e “o 

NG 

= max | Sado) e | (36) 
“o =o Zo 

e 
e 

ea 


ech respect to G. 


C. EXTENSION TO REDUCED-ORDER OBSERVERS 

So far, a new method to construct a full-order observers 
has been discussed. The method can easily be applied to 
design reduced-order observers also. Consider again the 


lWimear time-invariant system 


x(t) = A x(t) + B u(t) (37) 


| 
IQ 


y) xD (38) 


The reduced-order observer equation for this system is 


alt) = (Ao! = G Ap" 206) + (Ap, - CA)Z') G y(t) 


+ (Ap! = G A) y(t) + (Bs! wc B+') u(t) 


(39) 


20 





with 


w(t) = z(t) + G y(t) (40) 


where w(t) is the estimate of w(t). 

The characteristics of the observer system (39) depend 
on the coefficient matrix A,,' - G Ajo! which depends on the 
gan matrix G. Again to determine the gain matrix G, 


N 
consider the state error, w(t) - w(t). The observer error 


eau) is 

elt) 4 w(t) - w(t) (41) 
Differentiating equation (41) gives 

; ` : 

e(t) = w(t) - w(t) 


A ; 
Supo titutineg w(t) and w(t) and solving for e(t) yields 


/ 
E 2 (Aso iù Si Ay 5 ) e(t) (42) 
Letting F ê Ana! = G A,,', then (42) can be written as 
e(t) = Felt) (43) 


The solution to this equation is 
e(t) = et * e(o) (44) 


Notice that equation (44) has exactly the same form as 
Equation (27) except for the definition of the matrix F. 
ner sche results obtained for fuli-order observers 


apply to reduced-order observers as well. 


aa 





» 


DR SIMPLIFICATION FOR THE INFINITE-TIME CASE 

It will now be assumed that the upper limit on the 
Mco rmance messure integral is infinity. It will be shown 
Maat if Ca the problem of calculating v 8) can be 
simplified to the problem of solving a set of n(n+1)/2 
lbmeor algebrale equations. Let us state an important 
lemma and two theorems related to this simplification. 
ai NI SE Sep is to write an explicit solution for V(t). 


Lemma 1: The solution of the matrix differential 


equation 
V(t) = F v(t) + v(t) F+ Q (45) 


is 


E 
V(t) = g(t) V(O) dit) + I HE - +) Qt (t - +) dr (46) 
O 


where è(t) is the transition matrix for the system 
2 m 
z(t) = FE z(t) (47) 


Esso: Suppose that (46) is the solution of (45), and is 


rewritten as 


T 5 T T 
me) = p(t) V(O) o (t) + g(t) E ; g(-t) Q 4 (-t) at] 9 Ct) 


jit 


g(t) VCO) 9 (4) + o(t) L(t) o Ct) (48) 


where L(t) is defined as 


t 
A 9-1) Q 9 (1) at (49) 


L(t) Ê J 
0 


22 





Differentiating (48) gives 


W(t) = $lt) VO) 9 tt) + o(t) VO) olt) + olt) LCt) OTA) 


+ o(t) L(t) b(t) + lt) Totes) @ 6-01 62 (4) 
(50) 


The last term of equation (50) becomes Q, because 


POUR Oo eo (-t)] 0 (E) = 6(0) 96 (0) Tang 


Also 
p(t) = FÎ ọ(t) (51) 


merefore (50) can be written as 
V(t) = PT ot) V(0) di (t) + AML) V(O) $ (4) F 


+ pT o(t) L(t) T(t) RCE EE) We Heno 
(52) 


But from equation (18) 





a(t) VCO) T(t) = Vt) - o(t) Lit) f(t) (55) 
ubsvitutinae (53) into (52) gives 
W(t) = FY (v(t) - 86) L(t) g7(t)} + Cv(t) - $(t) L(t) 9 (6)] E 


t pi g(t) Ltt) ot (t) a) oT (t) F+Q (54) 


ran Tv ne equation (54) yields the final desired result 


v(t) = Fi v(t) +VF+Q 





The following well-known theorem gives a necessary and 
PeuilLetent condition for the error system (26) to be asymp- 
boulcally stable, that is, 


lim e(t) = 0 for all e(0). 


t>% 
Theorem 1: The linear time-invariant error system 
Sb) RE ct) (50) 


iS asymptotically stable if and only if, given any positive 
Simio Symmetric matrix @, there exists a positive 


Mete symmetric matrix K such that 


FK+KF=-0Q (56) 


ran 


Theorem 1 has been proved in [2]. This theorem can be 
used to determine whether or not the error system (26) is 


deniptovically stable for any choice of the gain matrix €. 





Wi-@icillowving taecrem applies speciiically to the case 


where têdê 


Theorem 2: Assume that the linear time-invariant system 
e(t) = F e(t) Op 


imemasympvOotically stable, and has a performance measure 


defined by 


t 
Ja f eltT)Qe(lr) dt (58) 


2l 








which is equivalent to 


J = (G) e (59) 


e V 
= — o 


ere V(t) satisfies 


V(t) = BÎ v(t) + V(t) F+Q 
with 
Y(0) = 0 
then 
v(t) = Fi v(t) + V(t) F+Q (60) 


WS a constant steady-state solution NERE ESCHE 


lim V(t) = Y, (a constant) 


t- 


Proof: By Lemma 1 with V(0) = 0, the solution to (60) is 


t 


vt) = So bt - 1) Q 6 (t - 1) at (61) 
Û 


, L 
Substituting o(t) = el Pana oTt) = el È into (61) and 


Meee properties of the matrix exponential, gives 


Ah JE T 
V(t) RCA teeny! az 


|| 

m 

| 
os 


pEr eon oe -F pi (62) 


2 








But the system (57) is asymptotically stable, so according 
Bezimeorem | there exists a positive definite symmetric 


matrix LE such that 


gizî. na mn (63) 


LL rs 
E — 


[F FT Ft 


WAL AN Ra A ce =Car) 


Since e T and pi commute, this can be written as 


pt t T 


_ E t A Aa -Ft Pt 
RECT) E ; dn le pe Te 
T T L 
e Ga î n 
—0 
0 
T ft T 
enni nee Ft 
ME E we e E E IV le 
T 
Ê Ft Ft 
RE me e (65) 
ie system (57) is stable, so 
Ft | rit 
lim e> = 10 ad lim e> = 0 (66) 
->00 ` r L>00 
Memce, substituting (66) into (65) 
lim V(t) = V (a constant) 


{t -+œ = 


iene Oves the theorem. 
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According to the result of Theorem 2, V(t) becomes zero 


ERR o Therefore, We is the solution of 


FFV+VF=-Q (67) 


which is obtained by setting V(t) = in equation 2) 
Equation (67) gives a set of n(n+1)/2 linear algebraic 
pellictdl'ons that can be solved for V after any positive definite 
Q has been selected for the performance measure (28). Hence, 
the problem of calculating LO becomes simple. The positive 
pera iteness ol the solution V can be checked by using 

Me ter's criterion [7] and if and only if V is positive 


definite the system is asymptotically stable. 


e 





IV SOLUTION OF THE OBSERVER DESIGN PROBLEM 


As mentioned before, the error equation needed to 
= 7 mine the gain matrix G has the same form for both the 
full-order and reduced-order observers. So one design 


method can be applied for both cases. 


A. DEVELOPMENT OF THE METHOD 


Ride arca that the gain matrix G be found which 





minimizes 
; ê CO m 
Jy = max | lame (To) 2 Gr at | 
2 0 
0) 
ca 
= max [ e Û WA ame (68) 
= <=" =o 7 
e 
Zo 
sce 
Letting min Jy = J,*, then 
J.* = min max [e î V(G) el (69) 
IL G o ¬0 Zo 
= e 
== 
lee l 
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Wii iS a minimax problem, and equation (69) is a real 
Suadratice form in Che There is a useful thecrem about the 
extremal properties of the eigenvalues of a real quadratic 
vm [1]. 

Theorem 3: The global maximum oí a real quadratic form 
on the unit hypersphere is equal to the largest eigenvalue 
of the quadratic form, and moreover the corresponding eigen- 
Repor 16 the vector drawn from the origin to the point on 


bil-Whyvpersphere where the quadratic form achieves its 


maximum. 
Therefore 
max BI by (CO) cee) (70) 
O a D 
e 
==) 
ee 


where 1, (6) is the largest eigenvalue of ve). Hence, now 
the problem is to minimize A,(G) with respect to the elements 


of G; that is, to-find the observer gain matrix G such that 


SE min Ay (6) Cie) 


To ensure that the observer is stable and does not have an 
unlimited bandwidth, the gain matrix G will be selected from 


the set 4 dei neg as 


gets 


1. All eigenvalues of F have negative real parts. 


A uno cai cenvalves ol = trace DL 


co 





where D is a specified positive constant. For a given value 
BRR tune largest eigenvalue is easily found by using a 
digital computer. The program required is simple because 
Ey efficient methods are available for finding the largest 
eigenvalue of a real symmetric matrix, and many computer 
facilities have several standard subroutines that can be 


ma for this purpose. 


B. ALGORITHM FOR COMPUTING THE GAIN MATRIX 

The steps in the computational procedure for finding the 
optimal gain matrix are given below. 

Guess G subject to the following requirements: 

Is stable. This can be done by assigning arbitrary 
eigenvalues with negative real parts to F and using Luenberger's 
method [4], for example. 

Bu ee > =D 


Zee 4 Minimization subroutine to minimize 





= x = 
ie ax 4 E e, A 


ee 


memvally, a penalty function is added to ensure that the 


trace F > -D constraint is satisfied, e.g., 


OF It gie lay 
J, + 


os piece aa DJS, ace 


tal 

lv 
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e, 
il 
ES 
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where a is a positive weighting factor. 
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ia computing: values of Já For che aio subroutine 
the following steps are carried out. 

aa solve pf ne = QO for = Vo 

b. Check to see that vr is positive de timite tt dit 
is continue, otherwise set ur = 10°, that 1s, a large 
meotuive number and return to the minimization subroutine. 


G If Vo TS positive definite it is the steady srate 


Solucion of 
: ta Te 
EG Sn NC) + Ve) Rae 


Next find the largest eigenvalue of LAN using a 
suitable subroutine. 
Gee Calculate the trace of F and compute the penalty 
Re & ON Y, 
ET ur = J, + Jo andererurmelostchesmininezatton 


al 


subroutine. 


In applying this algorithm to the examples given in this 
thesis the minimizations were accomplished using a subroutine 
"DIRECT" (NPS Computer Facility) that performs the pattern 
search method of Hooke and Jeeves, and the largest eigen- 
values were found using a subroutine "GIVHO" (NPS Computer 


Facility) that utilizes the Givens-Householder method. 
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RR A FULL-ORDER OBSERVER 

e uno se o thais eramplecis co illustrate the appli- 
selon oT the described algorithm to the design of a full- 
Rr Observer for a second-order system. 
a -—>—— Lo with the state- 


se +38 +2 





Consider the plant 


N 


variable representation 





x(t) = A x(t) + B u(t) 
0 il x (t) 0 
Ê + u(t) Cre) 
-2 -3 x(t) 1 
y(t) = C x(t) 
x (t) 
- [ı1 0] (139 
x (t) 


A full-order observer for this system is determined by 


specifying the observer gain matrix 


5] 


ES 
Il 


So 
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The observer system coefficient matrix is 


8¬ T 
F=A-GC= (74) 
-2-87 -3 


which has the characteristic equation 
A + (3+g,) DA + 38, + + 2=0 (75) 
61 1* 82 


FERE Crine the observer eigenvalues arbitrarily at -3.£j53 
gives the characteristic equation 


2 


ee) (X1+3-]3) = à" + 6A 4°18 = 0 (76) 


Equating the coefficients of like powers of 1 in (75) ana 


WHO) yields 


ja 
ji 


Mims the starting guessed value of the gain matrix G. If 
one arbitrarily chooses Q = I, and solves (67) for V the 


result is 


2 
UEA Ee ORO 
E 1 2 2 
kan 2 
12 + 228, + 08, + 28,85 + 6g, 
ê 3 - “8, - 81B> 
= AS (77) 
12 + 228, + 6g, + 28,85 + 6g, 
set See ee eo 
PE 


LD 2287 + 08, + 28185 + 68° 
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Due to the symmetry of V only three equations are needed. 
Ir the sum of the observer eigenvalues is to be limited to 
-10.0, for example, then the proposed algorithm gives the 
solution for the observer gain matrix as 

1220 


G= (78) 
-0.64063 


which gives the observer eigenvalues as 


-6.625 


1 


(79) 
2:315 
Eher selections of the weighting matrix Q and different 
lessons values for the trace of F the values obtained for 
Em mobserver gain matrix G and the corresponding eigenvalues 
are tabulated in Table 1. Looking at the table it can be 
seen that the two eigenvalues become farther apart from each 
other when the limiting value of the sum of the eigenvalues 
is decreased (assuming a fixed Q). 

To compare this full-order observer, which will be 
mMeterred to as the optimal observer, with other observers 
having different observer eigenvalues, consider two observers 
Moving eigenvalues selected as 

-5.0 -8.0 


Sg > NGE (80) 
-5.0 -2.0 


respectively. The corresponding observer gain matrices are 
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O ê 
Cala > G = (81) 
2.0 -7.0 
wespectively. 

For simulation purposes, the input u(t) was selected to 
Mema Unit step function. Assuming that the initial error 
states lie on the unit circle, the actual and the estimated 
states computed by using the optimal gains of equation (78) 
ape ea as functions of time on the graphs in Fig. 5 for 
rêv best initial error condition and in Fig. 6 for the worst 
iria! error condition. The estimate errors for the 
observers with three different observer eigenvalues are 
mented On one graph in Figs. 7 and 8 for the best and worst 
initial error conditions, respectively. The values of the 
performance measure for the observer error systems with 
different observer eigenvalues (different observer gains) 


are tabulated in Table 2 For che worsc inata! conditions, 









Worst. Initial Values of Perform- 
Conditions ance Measure 


Observer Gains. 
(Eigenvalues) 





= 0 0.16 
e(0) = 0.164 

| = -0.64063 ào = -0. 987 
oe! 
05505 


TABLE 2. PERFORMANCE VALUES FOR VARIOUS OBSERVER 
EIGENVALUES 
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B. A REDUCED-ORDER OBSERVER 

The purpose of this example is to illustrate the appli- 
cation of the algorithm to the design of a reduced-order 
ebserver for a third-order plant. 


Consider the plant characterized by the transfer function 


< 
^~ 


S) _ 1 
SI + 18º + 5942 


È 


WAN MD Lane can be represented by the state and output 


equations 
EO) rîn fito ai x, (6) 0 
EC) = 0 0 1 x(t) 4 Oi | Ge) Qo) 
x(t) -2 -5 ll BAV Û 
x(t) 
gd = [1 0 01 | 3,60) (83) 
x3(6) 


This plant is a third-order system with a single output so 
ER econd-order observer with arbitrary eigenvalues can be 
constructed. The output matrix C already has the same form 
as the one developed in the reduced-order observer theory. 
Hence, it is not necessary to change the variables, that SA 


for this case the matrices 
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The matrices A and B can be partitioned as 


0 1 1 0 0 
== = kas ma a am am — — = 
A =| 0 , 0 1 B=| 0 
i 
E Lins: i 


Therefore, the: reduced-order observer system matrix is 


281 E! 
Fz A GAGE (84) 
BE = 
62 


Der has the corresponding characteristic equation 


né + (4 + 8,)) 1 187 VE O (85) 


Suppose one selects the starting observer eigenvalues at 


A = -3, -2. This gives the characteristic equation 


ei) ree (86) 


Equating the coefficients of the same powers of A in (85) 


and (86) yields 


T0 


ES 
Ji 


-3.0 


This is the initial guessed value of the observer gain matrix 
ER TO the welgnting matrix Q is chosen arbitrarily as Q = ibe 


and equation (67) is solved for V the result is 
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36 + 4g, + lig, + A 





oz 2 
HO + 128, + n + 28,85 + 387 
y o L = 287 a 5780 
sl; 8g. 2 er 
E 
Gl es 2 
20 1 2 1 
22 


HO + 428, + 98» + 28,8, + 8g,“ 


RREO the Symmetry of V only three equations are needed. 
proposed algorithm gives the solution for the observer 


gain matrix as 


eV 
G = (88) 
-1.73438 
mame sum of the eigenvalues of F is limited to -30.0. 
ais gain matrix G gives the observer eigenvalues 
-25.851 
A = (89) 
-4.149 


o compare this optimally designed reduced-order 
observer with others having different observer eigenvalues, 


consider two observers with the eigenvalues 








-15.0 -20.0 
A = , A = (900 
-15.0 -10.0 


meepectively. The corresponding observer gain matrices are 


| 26.0 2630 
ce = 
116.0 91.0 


| 
li 


(ga) 


respectively. In Fig. 9 the actual and estimated states 
Meme the optimal gains are plotted for the best initial 
wre Conditions and in Fig. 10 for the worst initial error 
Hin LOS. It is assumed that the initial error condition 
ES on the unit circle. The estimate errors for the three 
Seeervers are plotted on one graph in Figs. 11 and 12 Tor 
mee best and worst initial error conditions respectively. 
Table 3 shows the values of the performance measure for 
observers with different eigenvalues for the worst initial 


cenditions. 


ER OBSERVATIONS 

The examples indicate that the observers designed Dy 
Being the proposed technique avoid large errors at a cost of 
having relatively small errors that slowly approach zero. 
ic is to be anticipated from the form of the performance 
Measure selected. 

Wo de sis approach is, however, Conservative im the 


sense that the resulting observers function effectively for 
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Values of Perform- 
ance Measure 


Worst Initial 
Conditions 





Observer Gains 
(Eigenvalues) 





TABLE 3. PERFORMANCE VALUES FOR DIFFERENT OBSERVER EIGENVALUES 


è 





selon in Which che initial conditions are worst. In 
PENNO OO enS  characteristic, performance for less 


dica lo initial conditions is compromised. 
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V. CONCLUSIONS 


In this thesis a new method for designing observers nas 
been developed. The optimal observer gains are determined 
by Solving a minimax problem. ‘The proposed algorithm is 
easily suited to digital computation and especially appli- 
fable to higher-order systems. 

The observers computed by using the algorithm have a 
PT Tdency to avoid large errors at the expense of small 
Dos mea tony decay. AS anticipated, the observers 
designed by the proposed technique are efficient for the 


Rest initial error conditions. 
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Revo Che initial state by 
e(tp) = o(t,) e(0) (A.5) 
Differentiating (A.4) with respect to t gives 


ot, nee or ro) 


Berl u. tana with the boundary condition 


f 


The performance measure for (A.1) is given by 


Es 


I=f Let) Gel) Ja WW 
0 


pepstitubing for e(t) using (A.2) gives 


Dr 


Jes el(0) #°(1) Q g(t) e(0) J at (A.8) 


Rearranging (A.8), one obtains 


t 


r 
JI = (0) LS HT) Q $lr) är Je(o) (A.9) 
0 
De fining 
E 
vit) As et(« - t) Q e(«- t) ar (4.10) 
then t 
M(0) = 4 Tê CA air (A.11) 
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Therefore, we have the final result (29) 


J = e!(0) V(0) e(0) (29) 


Differentiating (A.10) with respect to t by using Leibnitz's 


pule 
e t ê , 
Mit) = / et) Jar 
t 
T 
= $ (t - 1) Q t) 
È | | 
E [ o (1 - t) Q@ o(x - t) + oi (rt - t) QA dlı - t) Jr 


- Q (A.12) 


peostituting (A.6) for NG. - t) gives 

E Es 

W(t) = J | E se ea t)] 
t 


da Lo (1 - t) Q o(t - t) FI} dt - Q (A.13) 


Rearranging (A.13) gives 


ct 


; f m 
V(t) = - F i) @ (t= t) Q g(t — t) dr] 
IE 


“Lys or -t) Q (1 -t) ar] F- Q (A.14) 
t 


54 





pupetituting for Vit) using (A.10) gives equation (30) 


v(t) = - FE v(t) - v(t) F-Q (30) 


| 


Letting t = Es Mona aiyes equation (31) 


V(t.) (31) 


Il 
|o 


Taking the transpose of (A.10), the right side does not 


change, hence V(t) is symmetric. 
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